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ESTIMATION OF POND EVAPORATION 
IN THE GEORGIA COASTAL PLAIN 

By Richard H. McCuen^ and Loris E. Asmussen^ 

ABSTRACT 

Various techniques to estimate pond evaporation are compared on the basis of 
their ability to reproduce observed evaporation rates from a floating pan. Methods 
that use only one meteorological factor provided correlation coefficients that were 
only slightly smaller than those obtained with methods requiring values for two to 
four meteorological factors. Factor analysis indicated that radiation, humidity, 
and air temperature have similar effects on the variation in floating-pan evapora- 
tion; thus, it is not necessary to use all three factors for predicting evaporation 
rates. Correlation analysis, sensitivity analysis, and factor analysis indicated that 
variation in air speed has little effect on variation in evaporation. 


INTRODUCTION 

The mean annual rainfall for the 48-year period 
1923-70 recorded at the Coastal Plain Experiment 
Station, Tifton, Ga., was 47.24 inches. This was well 
distributed, with the months of October and 
November showing the lowest average rainfall of 
1.93 inches and July the highest of 6.03 inches. 
However, mean monthly rainfall is misleading be- 
cause there are several drought periods during the 
crop-growing season each year. Storage of suffi- 
cient water to meet irrigation requirements is the 
obvious solution; hence, estimates of pond evapora- 
tion in the Georgia Coastal Plain are necessary for 
estimation of pond storage requirements. 

The net evaporation rate is controlled by the 
meteorological conditions of the overlying air and 
the energy state of the evaporating surface layer. 
The quantity of thermal energy stored in a pond or 
lake can also be a controlling factor during certain 
periods of the year. Evaporation rates from open 
water surfaces have been estimated from mea- 
surements of water losses from pans (16)^ and 
sunken tanks (i), by using models developed from 


^Associate professor, Department of Civil Engineering, Uni- 
versity of Maryland, College Park, Md.; formerly postdoctoral 
fellow, Georgia Institute of Technology , at the time research was 
performed. 

^Geologist, Southeast Watershed Research Unit, Agricultural 
Research Service, U.S. Department of Agriculture, Tifton, Ga. 
31794 . 
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Cited” at the end of this publication. 


energy budgets (IS), and by empirical analysis of 
data (9). Methods that require measurements of 
water surface properties are rarely used because 
water surface data are difficult and expensive to 
collect. Thus, most evaporation estimation proce- 
dures are based on data that can be collected above 
the water surface. 

Hydrologic models are potentially useful for pre- 
dicting the output from a hydrologic process and for 
identifying the relative importance of variables that 
influence the hydrologic process being modeled. If 
the only objective of modeling is prediction, it may 
be reasonable to structure a model on only those 
variables that contribute statistically significant, 
physically meaningful information. Numerous 
methods involving meteorological factor * 

tor variables have been prop"^'"'"'^ 
evaporation rates. Howe^^'"' 
have nr- ^ 

r' 

examine the marginal value of four meteorological 
factors (radiation, air speed, relative humidity, and 
air temperature) for predicting evaporation rates in 
the Georgia Coastal Plain. 

PROCEDURES FOR 
ESTIMATING EVAPORATION 

Numerous methods have been proposed for es- 
timating evaporation. The prediction methods ex- 
amined herein are classified according to the data 
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required to estimate evaporation and whether or 
not the method requires empirical calibration of 
structural coefficients. The following three 
methods, none of which requires the estimation of 
structural coefficients, were examined: the U.S. 
Weather Bureau^ method (WB), the Penman 
method, and a modified fractional-evaporation 
equivalent method (MFEE). Additionally, the 
fractional- evaporation equivalent method (FEE) 
and four Mnear bivariate regression equations were 
calibrated fi^om observed data and then were used 
to predict daily evaporation. 

Weather Bureau Method 

The WB method of estimating evaporation was 
developed by Kohler, Nordenson, and Fox (7) fi"om 
an extensive analysis of data. The WB method re- 
quires observations of air temperature (T), radia- 
tion (i2), humidity (H), and air speed (XO. The 
method was originally presented in nomograph 
form; Lamoreux {8) provided the following func- 
tional representation: 

g[(T-2 12X0. 1024- 0.1066 In^J?)] _ Q QOQJ 

+ 0.025 (e, - eJ«-88(0.37 +0.004iy) 

X [0.025+ (r+398.36)"2 (4.7988x 10^«) €t-7482.6/(7'+398.36)]j -i(i) 

Here, and are the vapor pressures at dew point 
temperature and air temperature, respectively. 

Penman Model 

In a simplified form, an energy balance equates 
the net radiant-heat flux density (i2„) to the sum of 
the latent heat of evaporation {E) and the flux 
density of sensible heat {S), 

R„=E +S. (2) 

From aerodynamic considerations of turbulent 
transfer, the latent heat flux density is given by a 
form of Dalton’s law, 

^ =/o 00 (€„ - Ba), (3) 

and the sensible heat transfer is given by 

S = y/o 00 (To - T,), (4) 

where y is the psychrometric constant, /„ (V) is a 
wind function, e is the vapor pressure, T is the 
temperature, and the subscripts o and a refer to 
water surface and air quantities, respectively. 
Equations 2, 3, and 4 are the basis for the Penman 
model (IS, H.), given by 

E„=(^■R^ + y■ E,,) / (A + y), (5) 

where E„ is the evaporation estimate and A is the 
slope of the saturation vapor pressure curve at air 
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temperature (Ta). E^o is an empirical representa- 
tion of the latent heat flux density given by 

E,, = 0.36 (0.6 + WlOO) (e„ - ««). (6) 

where V is the air speed in miles per day. Stanhill et 
al. (18) showed that the ratio of net radiation to 
incoming solar radiation was approximately 58 per- 
cent for a small pond in Israel at approximately the 
same latitude as Tifton, Ga. Monteith and Szeicz 
(12) calculated a net radiation ratio of 0.53 for an 
open water surface. Davies (J^) indicated that net 
radiation (i?„) is, on the average, 55 percent of the 
measured radiation (R). Thus, Eg in equation 5 can 
be approximated by 

E, = [0.66 • A • R -y 

• 0. 35 (0.6 + y/100)(e„ - ej] /(A + y). (7) 

Modified Fractional-Evaporation 
Equivalent Method 

The MFEE method given by Linacre (10) relates 
the fraction of insolation (EIR) used for evaporation 
with the mean daily ambient temperature T^, in 
degrees centigrade, and the fraction relative hu- 
midity H by the following empirical nonlinear 
equation: 

EIR = 0.02 H (8) 

Linacre indicated that the coefficient 0.02 agrees 
favorably with values derived empirically for evap- 
oration from grass surfaces. Although the coeffi- 
cient for water surfaces would probably be smaller, 
a value for water was not available, and the coeffi- 
cient 0.02 is used herein. 

Fractional-Evaporation Equivalent 
Method 

The FEE method relates the fraction of insolation 
used for evaporation to air temperature (T^) by the 
linear equation 

EIR=g+fTa, ( 9 ) 

where E and R are evaporation and radiation, re- 
spectively, and/ and g are regression coefficients. 
The FEE method has, been used by Stephens and 
Stewart (19) in southern Florida and Lane (9) in the 
Western United States. Values of the regression 
coefficients were determined empirically for this 
study from the analysis of radiation, temperature, 
and pan evaporation data. 

Linear Regression With Meteorological 
Factors 

The following linear regression equation was 
used to relate individually each of the four 
meteorological factors (radiation, air speed, air 
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temperature, and relative humidity) to daily pan 
evaporation rate E : 

E = a + bF, (10) 

where F is the meteorological factor and a and b are 
the regression coefficients. The regression coeffi- 
cients were derived from the observed evaporation 
fi-om a floating pan as the criterion variable. 

RESULTS OF DATA ANALYSIS 

The Southeast Watershed Research Unit (Ag- 
ricultural Research Service) is conducting hydro- 
logic studies at a pond on the J. B. Walker farm, 5 
miles southwest of Tifton, Ga. Walker Pond has a 
maximum storage of 9.81 acre-feet and a maximum 
surface area of 2.56 acres. Daily evaporation is 
measured from a pan floating on the pond. The 
floating pan is partially shielded by a group of trees. 
The average tree height in the woodland (pine) east 
of the evaporation pan is 25 feet. There are about 
650 trees per acre in this plot. 

Meteorological data, including relative humidity, 
air temperature, air speed, and incoming radiation, 
are measured at the U.S. Weather Service office in 
Tifton, Ga. During fall and winter, the prevailing 
winds are fi"om the northwest and northeast. Dur- 
ing the remainder of the year, the prevailing winds 
are from the west and southwest. During a 
33-month period, 392 daily measurements of evap- 
oration from the floating pan and the four 
meteorological factors from the Weather Service 
office were made. 

Estimates of daily evaporation were derived by 
each of the previously described methods from the 
392 daily measurements of the 4 meteorological fac- 
tors and floating-pan evaporation data. The com- 
puted evaporation rates were compared with the 
corresponding measured values. The Pearson 
product-moment correlation coefficient R was used 
as a measure of the strength of the linear association 
between the observed and computed evaporation 
rates. The coefficient of determination is the 
proportion of the variance in the observed evapora- 
tion rates that can be attributed to its association 
with the computed rates. The resulting values ofR 
andiJ^ are given in table 1. 

The bivariate linear regression equations involv- 
ing humidity, radiation, air temperature, and air 
speed explained 38.5, 32.2, 31.2 and 1.5 percent of 
the total variation in the observed daily evapora- 
tion, respectively. The inadequacy of the air speed 
measurements for estimating evaporation is dem- 
onstrated by the negative value of the correlation 


Table 1. — Comparison of computed and observed 
evaporation rates 


Procedure R 

Weather Bureau 0.650 0.423 

Penman 626 .393 

Modified fractional-evaporation equivalent 623 .388 

Fractional-evaporation equivalent 579 .334 

Linear regression: Predictor variable: 

Humidity -.620 .385 

Radiation 567 .322 

Temperature 559 .312 

Airspeed ..-.122 .015 


coefficient. A negative value of R for such a rela- 
tionship indicates that increases in air speed are 
accompanied by decreases in evaporation. Such a 
result would not be expected. 

The FEE method, which uses two meteorological 
factors (radiation and air temperature) to estimate 
evaporation, explained 33.4 percent of the total 
variation, only 1.2 and 2.2 percent greater than 
the percentage of variation explained by the bi- 
variate regression equations for radiation and air 
temperature. This suggests that radiation and air 
temperature are intercorrelated and tends to ex- 
plain similar variation in evaporation rates. 

The MFEE method, equation 8, explained 38.8 
percent of the total variation in the observed evap- 
oration rates, which is only 0.3 percent greater than 
the value resulting from the bivariate regression 
equation for humidity. However, the coefficients 
for the bivariate regression equation involving 
humidity were derived from data collected at Tif- 
ton, -Ga., while the coefficients for the MFEE 
method were firom the developments of a previous 
investigator (id). Thus, the difference of 0.3 percent 
in explained variation is not necessarily an indica- 
tion that humidity, air temperature, and radiation 
explain similar variation in evaporation. Also, the 
coefficients of the MFEE method are possibly in- 
adequate for estimating evaporation from ponds in 
the southern Coastal Plain. 

The WB and Penman methods, which require 
values of four meteorological factors, explained 42.3 
and 39.3 percent of the total variation, respectively. 
The MFEE method, which involves only three 
meteorological variables and does not require cali- 
bration of structural coefficients, explained 38.8 
percent of the variation. Thus, air speed, which is 
used with the WB and Penman methods but not 
the MFEE method, does not significantly increase 
the percentage of explained variation. Such a re- 
sult was expected, since the correlation coefficient 
for the bivariate regression equation involving 
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‘ ■: .i:.- ".cry .-rr.ai; and not statistically 

.-r coefficients for the 

n invordng humidity 
fr ::: ' f data, this equation ex- 

'c!ven: less variation than 

’ d-r.r:Uv:: metir^ds. respectively. Thus, 

' ' .c • cxar:::!:e the possibility that 

‘ ’■ f ra hat:-^!:, ah* tenq^erature, and 

t: j.d:: <“n:;hir variation in evaporation 
‘ tnat the coiiectiun of data for 
* ‘ ' V ‘ “ V . - 1 c !“ . vica. fact ors !*epresent s a dupli- 
**' ^ r!nat.‘„n ar/.: effort for prediction of 

' V rates: a ractcr analysis will be used 
**vaa..:.o tr> interir nv.ations of the meteoro- 
t./ . r: a .ater section. 

SENSITIVITY ANALYSIS 

^ '* ' '^‘tridents are used to measure the 

; t.'-.'-'-ciatrsn between variables. The Pear- 
, i ’■ :o''t-:rr men: correlation coefficient is de- 


- Ah ,i = F™:}h and /I, is the number of. lyi/ 
of the correlation coefficient rep- 
;:**"i^* of the variance in a criterion 

' ■ tx - tr ibut ed to its linear reeu'ession on a 

i rear: !wV. 

^ analyse is used to define the relation- 

i: variables. There is a strong sti’uctural 
.'miLarity between eon-elation analysis and regi-es- 
M‘ .n analysis i / r i. The linear regression coefficient b 
can he computed from the correlation coefficient 
ami The standai-d deviation of Z and Y: 

*=2 /lx,^=RSJS,. ( 12 ) 


relationship between linear regression and sensitiv- 
ity. Therefore, estimates of sensitivity can be used 
to determine a priori the factors that might account 
for the most variation in the dependent variables. 

The value of sensitivity computed by the deriva- 
tive approach of equation 13 is not invariant to the 
magnitude of either the dependent or independent 
variable. Relative sensitivity R,, is defined as the 
ratio of the percentage change in evaporation 
(ilEIE) to the percentage change in the magnitude 
of the meteorological factor (AF/F): 

R .= ^=(^)( l \. ( 14 ) 

IF/F \ aF /\ e / 

Since the relative sensitivity is a function of F 
and Z, it can be used for comparing the sensitivity 
of the different meteorological factors that in- 
fluence evaporation. A negative value of R, 
indicates that as the value of the meteorological 
factor increases evaporation will decrease. 

In order to relate, relative sensitivity values to the 
influence of the meteorological factors in estimating 
evaporation rates, it is necessary to specify a model 
of the evaporation process from w'hichiS, values can 
be derived. Use of a model based on the underlying 
ph^-sical principles w-ill provide /?, values that are 
indicative of the relative influence of variables in the 
evaporation process. Thus, the Penman model was 
selected because it is based on a simplified energy 
balance equation. With average monthly values (20) 
of meteorological factors and the Penman model 
monthly values of relative sensitivity were deter- 
mmed (table 2). 

T^e mean annual relative sensitivity values for 
the four meteorological factors indicate that humid- 

Table 2.— Monthly relative sensitivity values for 
Tifton, Ga. 


'^’ouation 12 suggests that for values ofie near 1, 
A will produce cnmnavativelv large 


Mean daily 
Month evaporation 
(inch/day) 

JanTTATabsi 

101 

132 

... .172 
... .196 
... .199 
... .184 
... .180 
... .148 

119 

092 

073 

Mean T ^ ~ 


Relative sensitivity 
of the meteorological factor 


Humidity Radiation 

-1.034W^ 0655 

-•811 .258 .652 

-.672 .224 . 698 

-.569 .199 .732 

-.542 .189 .745 

-.521 .158 .777 

-.641 .128 .819 

-•656 .131 . 815 ' 

-•763 .152 .7g5 

-•805 .189 .732 

-.846 .232 .671 

- 1.184 . 277 . 626 

-.753 .199 720 


Temperature 


0.661 

.688 

.637 

.621 

.627 

.589 

.484 

.495 

.554 

.629 

.677 

.615 



ty is the most sensitive factor and that radiation 
ind temperature are slightly less sensitive. The 
nean value for air speed indicates that a 
. 0-percent change in air speed produces a 2-percent 
;hange in evaporation. Linsley, Kohler, and 
^aulhus (11) indicated that a 10-percent change in 
lir speed causes a change in evaporation of 1 to 3 
)ercent. Thus, the relative sensitivity values de- 
lved herein appear reasonable. A ranking of the 
nean sensitivity of the four factors corresponds 
ixactly to the ranking of correlation coefficients for 
he bivariate regression equations of table 1. In 
>oth rankings, there is a very noticeable difference 
>etween air speed and the three remaining me- 
eorological factors. 

The monthly R, values of table 2 indicate that 
vaporation rates are most sensitive to variation in 
tumidity during winter. During summer, variation 
n radiation produces a greater change in evapora- 
ion than does variation in humidity. From April to 
une, both temperature and radiation are more sen- 
itive than humidity. Thus, the values of relative 
ensitivity are potentially useful for indicating the 
elative importance of the factors that control the 
vaporation process and the contribution of each 
actor in the reproduction of measured evaporation 
ates. 

FACTOR ANALYSIS OF PAN 
EVAPORATION AND 
METEOROLOGICAL DATA 

Factor analysis, a multivariate statistical tech- 
que, can be used to examine the intercorrelations 
stween the observed meteorological and pan 
vaporation data. Such an analysis is desirable be- 
Luse it can lead to a reduction in the number of 
ndependent'' variables required to predict a “de- 
mdent^’ variable. Elimination of predictor vari- 
)les that explain similar variation in the criterion 
n reduce the cost of future data collection. A re- 
iction in the number of predictor variables may 
duce prediction ability somewhat, but for vari- 
)les that are highly intercorrelated, the reduction 
often insignificant. 

Factor analysis transforms the correlation matrix 
the criterion and predictor variables into a set of 
iriates called factors. These factors are linear 
mbinations of the original variables and have the 
[vantage of being uncorrelated. The coefficients of 
e fiictors are called factor loadings. 

The correlation matrix of table 3 was derived 
3m 392 daily measurements of four meteorological 
ctors and evaporation from a floating pan. With 


the exception of associations involving air speed, 
the correlation coefficients between each pair of 
variables exceeded 0.5. The correlation coefficients 
between air speed and each of the other four vari- 
ables are between zero and -0.3. But the correla- 
tion matrix by itself does not provide the means for 
determining whether or not the four variates other 
than air speed explain similar information. 

A factor analysis of the correlation matrix of table 
3 produced the eigenvalues and factor pattern of 
table 4. In the first factor all of the factor loadings 
are large except the factor loading for air speed. The 
factor loading for air speed is large in the second 
factor while the other four factor loadings are small. 
In the remaining tliree factors, all the loadings are 
comparatively small. 

Table 3. — Correlation matrix 


Pan Air Radia- 

Variable evap. Temp. speed tion VPD 

(1) (2) (3) (4) (5) 

(1) Pan 

evaporation . .1.0000 0.5587 - 0.1224 0.5667 0.6203 

(2) Temperature... .5587 1.0000 -.2914 .5282 .8240 

(3) Air speed -.1224 -.2194 1.0000 -.1981 -.2939 

(4) Radiation 5667 .5282 -.1981 1.0000 .7237 

(5) Vapor pressure 

deficit 

(VPD) 6203 .8240 -.2939 .7237 1.0000 


Table 4. — Eigenvalues and corresponding 
eigenvectors 


Eigenvalues for eigenvector of — 

Variable "“ToT“^.93~”™~Td49 Cm 7m 

(1) (2) (3) (4) (5) 

(1) Pan evaporation .. . 0.779 0.224 -0.305 -0.491 0.002 

(2) Temperature 854 .057 .478 -.051 .190 

(3) Ail- speed -.358 .925 .089 .085 -.021 

(4) Radiation 817 .101 -.366 .420 .110 

(5) Vapor pressure 

deficit 937 .011 .171 .121 -.279 


Although methods exist for determining the 
number of statistically significant factors (2, S, 15)^ 
a universally accepted standard for when to stop 
factoring has not been developed (5). Furthermore, 
information, even though statistically significant, 
may have no practical value (5). Therefore, a 
method proposed by Kaiser (6) through extensive 
data analysis was used to determine the number of 
significant factors. Kaiser recommended that the 
factors corresponding to eigenvalues of the correla- 
tion matrix that are greater than 1 could be consid- 
ered significant. 

The eigen valuaiS' of table 4, determined by 
Kaiser's method, indicate that only the first factor is 
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significant. The large factor loadings in the first 
factor, which correspond to radiation, vapor pres- 
sure deficit, air temperature, and pan evaporation, 
suggest that observations of the three meteorologi- 
cal and pan evaporation variables contain similar 
information even though the four variables are dis- 
similar. If the three meteorological factors explain 
similar variation in pan evaporation, then a predic- 
tion equation including more than one of the three 
meteorological factors should not be expected to 
provide significantly greater correlation than a pre- 
diction equation involving only one of the variables. 
Comparison of the correlation coefficients of table 1 
confirms the findings of the factor analysis and indi- 
cates that the methods involving more than one 
meteorological factor provide correlation coeffi- 
cients only slightly higher than those obtained with 
the bivariate regression equations. Furthermore, 
the very low correlation between air speed and 
evaporation supports Kaiser’s empirical rule for de- 
termining the number of significant factors. The 
sensitivity analysis also indicated that variation in 
air speed had little effect on variation in evaporation 
rates. 

CONCLUSIONS 

Since true rates of evaporation from the pond 
surface were not available, the observed floating- 
pan evaporation rates are used to represent them. 
Such an assumption seemed reasonable since the 
mean annual evaporation from the floating pan was 
only about 9 percent less than the annual rate re- 
ported for the Tifton area by the U.S. Department 
of Commerce {20). During fall and winter, air moves 
primarily out of the northwest and northeast. The 
wooded area north of Walker Pond could signifi- 
cantly influence evaporation rates and thus account 
for, in part, the difference between the observed 
and computed mean annual evaporation. The mean 
observed May-to-October evaporation expressed as 
a percentage of the total annual evaporation dif- 
fered from the reported evaporation by only 2 
percent. From May to October the wind moves 
primarily from the west and southwest. During 
this period the wooded area on the northern bank 
of the pond would not significantly influence evap- 
oration rates. Furthermore, a recent study {1) in- 
dicated that the average ratio of evaporation from 
Lake Hefner to evaporation from tanks with diam- 
eters of 9 and 15 feet and located approximately 
100 feet from the lake is 1.01. 

As indicated previously, an understanding of the 
evaporation process can provide the means for ob- 
:aimng better estimates of evaporation rates. Such 


estimates of evaporation can then be used to pro- 
vide more accurate estimates of the storage re- 
quired for irrigation supplies during drought. A 
sensitivity analysis was used to examine the rela- 
tive importance in the Georgia Coastal Plain of 
four meteorological variables thought to influence 
evaporation. Mean annual relative sensitivity 
values indicated that humidity, radiation, and air 
temperature are of similar importance, while vari- 
ation in air speed has a much smaller influence on 
evaporation rates. However, mean monthly rela- 
tive sensitivity values indicated that humidity is 
the controlling factor in winter, while radiation 
exerts a greater influence on evaporation rates in 
summer. 

The correlation analysis of table 1 indicates that 
the WB method is potentially capable of providing 
the best estimates of evaporation. If a method for 
estimating net radiation is available, the Penman 
method also provides reasonable estimates of daily 
evaporation. However, the correlation coefficients 
obtained with three of the bivariate regression 
equations were only slightly less than the correla- 
tion coefficients obtained with the WB and Penman 
equations; this suggests that the use of more than 
one meteorological factor is redundant when the 
objective of modeling is prediction. 

Whereas the coefficients of the WB and Penman 
methods were derived at other locations, the re- 
gression coefficients for the bivariate regression 
equations were derived from data collected at Tif- 
ton, Ga. Therefore, a factor analysis- was used to 
examine the intercorrelations between the me- 
teorological factors and floating-pan evaporation 
rates. The factor analysis supported the conclusions 
of the correlation analysis. Specifically, the factor 
analysis indicated that radiation, air temperature, 
vapor pressure deficit, and pan evaporation explain 
similar variation, and thus, for the Georgia Coastal 
Plain, the use of more than one of the three me- 
teorological factors (vapor pressure deficit, radi- 
ation, and air temperature) may be redundant for 
predicting evaporation. 

The correlation analysis, sensitivity analysis, and 
factor analysis indicated that evaporation rates are 
only slightly influenced by variation in air speed. 
The sensitivity analysis showed that a variation in 
air speed of 10 percent results in a 2-percent change 
in evaporation. Thus, the very low correlation coef- 
ficient, — 0. 122, between air speed and pan evapora- 
tion should have been expected. The factor analysis 
also indicated that air speed was not a significant 
factor in estimating floating-pan evaporation. 



LITERATURE CITED 

(1) Allen, J. B., and Crow, F. R. 1971. Predicting lake evapora- 

tion by performance of evaporation ponds, pans and 
tanks. Trans. ASAE 14(3); 458-463. 

(2) Bartlett, M. S. 1950. Tests of significance in factor analysis. 

Br. J. Psychol., Stat. Sect. 3; 77-85. 

(3) 1951. A further note of tests of significance in factor 

analysis. Br. J. Psychol., Stat. Sect. 4: 1-2. 

(4) Davies, J. A. 1967. A note on the relationship between net 

radiation and solar radiation. Q. J. R. Meteorol. Soc. 93: 
109-116. 

(5) Harman, H. H. 1960. Modern factor analysis. 474 pp. Uni- 

versity of Chicago Press, Chicago. 

(6) Kaiser, H. F. 1960. The application of electronic computers 

to factor analysis. Educ. Psychol. Meas. 20(1): 141-151. 

(7) Kohler, M. H., Nordenson, T. J., and Fox, W. E. 1955. 

Evaporation from pans and lakes. U.S. Weather Bur. 
Res. Pap. 38, 53 pp. 

(8) Lamoreux, W. M. 1962. Modern evaporation formulas 

adapted to computer use. Mon. Weather Rev., vol. 90, 
January. 

(9) Lane, R. K. 1964. Estimating evaporation from insolation. 

J. Hydraul. Div., Proc. Am. Soc. Civ. Eng. 90(HY5): 
33-41. 

(10) Linacre, E. T. 1967. Climate and the evaporation from 
crops. J. Irrig. Drain. Div., Proc. Am. Soc. Civ. Eng. 93 
(IR4): 61-79. 


(11) Linsley, R. K., Kohler, M, A., and Paulhus, J. L. H. 1958. 

Hydrology for engineers. 340 pp. McGraw-Hill, New 
York. 

(12) Monteith, J. L. , and Szeicz, G. 1962. Radiative temperature 

in the heat balance of natural surfaces. Q. J. R. Meteoroli 
Soc. 88: 496-507. 

(13) Penman, H. L. 1948. Natural evaporation for open water, 

bare soil, and grass. Proc. R. Soc. London, Ser. A, 193: 
120-145. 

(14) 1956. Estimating evaporation. Trans. Am. Geophys. 

Union 37: 43-46. 

(15) Rao, C. R. 1955. Estimation and tests of significance in 

factor analysis. Psychology 20: 93-111. 

(16) Shulyakovskiy, L. G. 1969. Formula for computing evap- 

oration with allowance for the temperature of the free 
water surface. Sov. Hydrol. Sel. Pap. 6; 566-573. 

(17) Snyder, W. M. 1971. Summary and evaluation of methods 

for detecting hydrologic changes. J. Hydrol. 12: 311-338. 

(18) Stanhill, G., Hofstede, G. W., and Kalma, J. D. 1966. Radi- 

ation balance of natural and agricultural vegetation. Q. J . 
R. Meteor. Soc. 92: 128-140. 

(19) Stephens, J. C., and Stewart, E. H. 1963. A comparison of 

procedures for computing evaporation and evapo transpi- 
ration. Int. Assoc. Sci. Hydrol. Pub. 62, pp. 123-133. 

(20) U.S. Department of Commerce. 1968. Climatic atlas of the 

United States. 104 pp. Washington, D.C. 


*l796-G.P.0.-1750-s/67l -583/46 



